Ehrenfest time dependence of quantum transport corrections and spectral statistics 
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The Ehrenfest time scale in quantum transport separates essentially classical propagation from 
wave interference and here we consider its effect on the transmission and reflection through quantum 
dots. In particular we calculate the Ehrenfest time dependence of the next to leading order quantum 
corrections to the transmission and reflection for dc- and ac-transport and check that our results 
are consistent with current conservation relations. Looking as well at spectral statistics in closed 
systems, we finally demonstrate how the contributions analyzed here imply changes in the calculation 
given in [P. W. Brouwer, S. Rahav and C. Tian, Phys. Rev. E 74, 066208 (2006)] of the next to 
leading order of the spectral form factor. Our semiclassical result coincides with the result obtained 
in [C. Tian and A. I. Larkin, Phys. Rev. B 70, 035305 (2004)] by field-theoretical methods. 
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I. INTRODUCTION 

Chaotic quantum systems are expected [l[ to show uni- 
versal behavior that can be described by Random Matrix 
Theory (RMT) 0. After this was conjectured [H the 
challenge was to justify and dynamically understand the 
relation between chaotic systems and RMT. Here semi- 
classical methods have proved to be very successful 0-0| ■ 
These are based on asymptotic expansions of the quan- 
tum propagator, the Green function and its trace, which 
consist of sums over classical trajectories ||. Each sum 
contains, along with prefactors determined by the classi- 
cal dynamics, phases determined by the classical actions 
of the trajectories which allow for interference effects. 

On one front, these semiclassical methods were applied 
to study spectral properties of closed systems. Here one 
considers for example the spectral autocorrelation func- 
tion 

K{u) = 2nh ( Posc (E + huj/2) Posc (E - fiw/2)) , (1) 

defined as the energy averaged correlation function of the 
oscillating parts of two densities of states p osc (E) with 
an energy difference Hu>. Here and in the following (...) 
denotes an average over a classically small but quantum 
mechanically large energy window AE. Semiclassics now 
enters by replacing the spectral densities by their semi- 
classical expression in terms of a sum over periodic orbits 
given by the Gutzwiller trace formula Q 



(E) 



(i/h)S~,(E) 



(2) 



for h — > with A 1 the stability amplitudes (for their 
exact form see [8j for example) and S 7 (E) the classical 
actions of the periodic orbits and denoting the real 
part. Linearizing the actions around the energy E finally 
yields when defining 




FIG. I: Correlated orbits analyzed by Sieber and Richter [f| 
which differ in an encounter and lead to the first off-diagonal 
correction to the spectral form factor. 



for 



K(lu) ~ 2&K (u) 



(4) 



with T 1 the period of the orbit 7. Often the spectral 
form factor K(t), which is Fourier transform of K(u), is 
considered 



K(t) 
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(5) 
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(iu,/2)(T 1 (E)+T 1 ,(E)) 



(3) 



Importantly, the expressions for K(oS) and K(r) os- 
cillate rapidly depending on E; dominant contributions 
will thus result from trajectories with very similar ac- 
tions. For example, the diagonal contribution, i.e. 7 = 7' 
with equal action, was first studied in 1985 in (3), using 
the sum rule of [4(. It yielded the leading order RMT- 
prediction in 1/ (iuS) for the spectral autocorrelation func- 
tion and thus the leading order in r for the spectral form 
factor. Off-diagonal contributions were first taken into 
account by Sieber and Richter in 2001 [j| by consider- 
ing two orbits essentially differing from each other in 
a encounter region where the two orbits are differently 
connected, see Figure [TJ This work could later be ex- 
tended and formalized yielding the RMT results to arbi- 
trary high order in powers of r [6, 7] and also extended 
to other quantities characterizing the spectral properties 
of a system [§]. 

On another front, for open systems the conduc- 
tance was analyzed semiclassically within the Landauer- 
Biittiker approach 10] to transport. In particular we 
imagine a system connected to two leads carrying Ni and 
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N% channels respectively (with a total N = iVi + A^)- 
The conductance is related to the transmission matrix 
elements t a ,p via the Landauer-Biittiker approach and 
these matrix elements can be approximated semiclassi- 
cally (for an overview see [ll|) by 




B e (i/fi)(S^(E)) 



(6) 



with ta t p characterizing the transitions between the 
modes /3 and a, and 7 the classical scattering trajectories 
that connect those modes in the two different leads. The 
Heisenberg time is defined as Th = 2i:hd{E) with d(E) 
the mean spectral density of the considered closed sys- 
tem. For the exact form of the stability amplitudes £? 7 



and classical actions S 1 (E) see e.g. (llj. A corresponding 
expression in terms of trajectories connecting one lead to 
itself holds for reflection amplitudes r a ,p- Using eq. © 
we obtain for the transmission T characterizing the con- 
ductance 



T= (Tr(tt f )) 




The sum runs over all paths 7 and 7' which connect the 
two leads. A similar expression also holds for the re- 
flection R = (Tr (ttt ) ). To be precise we will use R 
to denote the reflection into lead 1 in the following, but 
the expression for the reflection into lead 2 just follows by 
swapping N% and A^. The first nondiagonal (i.e. next or- 
der in inverse channel number 1/AQ contribution for the 
conductance was analyzed in [12| yielding again a result 
consistent with RMT. The extension of the conductance 
to arbitrary high order was performed in [l3| and the 
authors later treated the shot noise [lij and other cor- 
relation functions like the conductance variance in [l5j . 
For this the authors built on their work on closed sys- 
tems Q and noticed in particular that the diagrams of 
correlated pairs of scattering trajectories that appear for 
the conductance can be created by cutting the pairs of 
periodic orbits that contribute to the spectral form fac- 
tor once and moving the cut ends to the leads. Likewise, 
the diagrams of trajectory quadruplets that appear for 
the conductance variance can be obtained by cutting the 
periodic orbit pairs exactly twice. 

Up to now we only discussed one application of 
these semiclassical techniques, the confirmation of RMT- 
results. However it is also possible to predict effects 
away from this arena, i.e. the behavior of chaotic sys- 
tems for a finite Ehrenfest time. The Ehrenfest time 
t~e = (1/A) \n(E / (XH\}\ more generally defined as a time 
proportional to In h [16[ ; is the time needed for a wave- 
packet to reach a size such that it can no longer be 
described as a single classical particle. The Ehrenfest 
time thus separates the free evolution of wave packets 
that follow essentially the classical dynamics from the 
evolution on larger time scales where wave interference 
becomes dominant. Including the Ehrenfest time into 




FIG. 2: Diagram occurring in the calculation of the reflec- 
tion covariance (or the conductance variance) containing two 
orbits surrounding a central periodic orbit. The fringes are 
marked by (black) vertical lines perpendicular to the trajec- 
tories. Partner orbits are not shown. 



the calculation of the quantities presented above started 
with the pioneering work [ItJ that calculated the first 
quantum correction to the energy-averaged transmission. 
This analysis was extended to reflection [l^, [l9[ includ- 
ing also a distinction between different Ehrenfest times 
[18j. An exponential suppression of this quantum cor- 
rection proportional to qt Te I Td was observed involving 
the dwell time td, the average time the particle stays 
inside an open billiard. Furthermore the Ehrenfest time 
dependent behavior of other transport quantities soon 
followed: the independence of (the leading order of) the 
universal conductance fluctuations was obtained in [l9j . 
the shot noise and Fano factor were found to be exponen- 
tially suppressed like the averaged transmission in [20| . 
and the behavior of a third order correlation function was 
derived in [2lj |. 

Of these it is the treatment of the conductance variance 
fl9j we are particularly interested in here. Because of the 
unitarity of the scattering matrix this is equal to the re- 
flection covariance, which turns out to be slightly simpler 
to treat semiclassically, and the authors found one im- 
portant contribution was given by a diagram like in Fig- 
ure [21 There two trajectories (one from either lead) ap- 
proach a trapped periodic orbit with one winding around 
it an extra time. Partner trajectories (not shown) can be 
found which follow those trajectories almost exactly but 
where one winding is exchanged between the two trajec- 
tories leading to a quadruplet of trajectories with a small 
action difference and a contribution in the semiclassical 
limit. Such a contribution vanishes when the Ehrenfest 
time goes to and can be seen to contain the discrete dia- 
gram types considered previously without Ehrenfest time 
15] (which then naturally sum to 0). Although this con- 
tribution vanishes, similar periodic orbit encounters can 
contribute in other situations 0, [22j when the Ehrenfest 
time is 0. 

Combined with another diagram which does not in- 
volve a periodic orbit encounter, [l9j showed the inde- 
pendence of the conductance variance of the Ehrenfest 
time. Later these techniques were applied to the spectral 



3 



autocorrelation function and the spectral form factor for 
closed systems in (23|. This allowed the authors of [23] to 
obtain the first quantum correction in the case with and 
without time reversal symmetry, but they found a dis- 
crepancy with the field-theoretical result [24[ obtained 
using effective RMT, a phenomenological approach to 
mimic the Ehrenfest time behavior in the RMT frame- 
work. We investigate this discrepancy here in this article 
and show how a hierarchy of diagram possibilities (e.g. re- 
versing the cutting of periodic orbits to create diagrams) 
restores the consistency in the semiclassical treatment. 

Currently all the Ehrenfest time approaches (described 
above) are restricted to very low order in the inverse 
channel number for the transmission and in 1/lo or r 
for the spectral autocorrelation function or form factor. 
A calculation of the corrections to infinite order, as has 
been performed in the case of vanishing Ehrenfest time 
te/td —> 0, is still lacking. We want in this paper to 
make a step towards filling this gap. More precisely we 
consider in the section II the next-to-leading order quan- 
tum correction to the transmission and reflection in the 
case of the dc-transport with and without time reversal 
symmetry. We then check the unitarity of our result, i.e. 
that T and R add up to a constant (Ni) at the consid- 
ered order. In section III we extend the results of section 
II to ac-transport and then check that corrections to the 
closely related Wigner time delay are indeed zero at the 
order considered. In section IV we apply our previous 
results to closed systems to obtain for the spectral form 
factor with Ehrenfest time a result consistent with the 
field-theoretical prediction [24| and finally conclude. 
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FIG. 3: Example of an orbit (and its partner shown dashed), 
considered in [l3j], that contributes to the transmission for 
systems without time reversal symmetry. A central periodic 
orbit (dashed-dotted) can be identified. 



off-diagonal corrections cancel. Of course as both involve 
a closed loop which is traversed in two different directions 
by the trajectory and its partner, they do not exist and 
can yield no contribution when time reversal symmetry 
is absent. 

But it is the next order contributions we are partic- 
ularly interested in, and we start with the simpler case 
where the scattering system does not have time reversal 
symmetry. 



II. TRANSMISSION AND REFLECTION 

Before we turn to the form factor later we remain with 
quantum transport and consider the transmission and re- 
flection. The leading order contribution (in inverse chan- 
nel number l/N) to the transmission and reflection re- 
sults from the diagonal approximation (pairing 7 = 7' in 
(O). The calculation of this contribution can be found 
for example in [25| and the result is independent of the 
Ehrenfest time and is of order N. The next order in in- 
verse channel number (i.e. of order 1) results from the pe- 
riodic orbit pairs shown in Figure [TJ where one of the two 
loops (that which is traversed in the same direction by 
both orbits) is cut open and the two ends are brought to 
the two openings [H, EH and its contribution is damped 
exponentially with the Ehrenfest time [TtJ • For the re- 
flection an additional possibility arises, called coherent 
backscattering, and which can be created by cutting the 
orbits in Figure [T] in half (keeping the half traversed in 
different directions by the orbits) and moving what is 
left of the encounter to the lead. As there is still the 
remnant of the encounter, this case is also suppressed 
exponentially with the Ehrenfest time [TH, [r| . This de- 
pendence is essential for the unitarity of the scattering so 
that if we sum the transmission and the reflection these 



A. No time reversal symmetry 

The l/N order contribution results from orbits with 
two encounters with itself [13| . see Fig.[3j We can see that 
there is a central periodic orbit through the two encoun- 
ters. This fact is essential for the Ehrenfest time depen- 
dence and simplifies treating the different cases. Depend- 
ing on how much these encounters overlap (i.e. depend- 
ing on the lengths of the links ti and £5 in Fig. [3]), one 
distinguishes in the case of no overlap two independent 
2-encounters (i.e. encounters involving 2 orbit stretches), 
in the case the two 2-encounters overlap at one of their 
ends (shrinking t± or £5 say) a 3-encounter and in the case 
the two 2-encounters overlap at both ends (shrinking t 4 
and ts) an encounter fully surrounding the periodic orbit. 
Although we mentioned up to now only one contained pe- 
riodic orbit shown dashed-dotted in Fig. |3l there are two 
in total: one built up by £4 and ts, the other by t% and £4. 
In the following calculations we choose either as we actu- 
ally treat this configuration as an orbit meeting a central 
periodic orbit twice, see Fig. |4j This procedure counts 
every configuration twice; this overcounting factor ac- 
counts for the fact that for fixed orbit parts, i.e. for fixed 
dashed-dotted periodic orbit and fixed orbit encountering 
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it in Fig. [4] we have two possibilities to construct an orbit 
pair: the original orbit can surround the dashed-dotted 
orbit once more than its partner during either the first or 
the second encounter. The two possibilities correspond 
to swapping the original orbit and its partner, and both 
terms are included in the sum over orbits in ([7]). The 
equivalence between chosing a different central periodic 
orbit and swapping 7 and 7' can be seen for example 
by following the orbits in Fig. [5] As we later sum over 
all possible central periodic orbits, we fix 7 as having one 
traversal fewer than 7' during its first encounter with the 
central periodic orbit and one traversal more during its 
second encounter. 

The two encounters of the orbit with the central 
periodic orbit have to be independent, because other- 
wise there exists no connected partner with a small but 
nonzero action difference. That means the orbit has to 
decorrelate from the central periodic orbit. Then this or- 
bit has to become ergodic before returning, therefore an 
encounter time, that is of the order of the Ehrenfest time 
is required, so the stretch away from the periodic orbit in 
Fig.[3]or the top loop in Fig.[4]must be of positive length. 
The total orbit has thus to be longer than the sum of the 
two durations of the encounters with the central periodic 
orbit, i onc ,i + icnc,2- This excludes the case that both £3 
and £5 in Fig. [3] get so short that both stretches do not 
decorrelate from the central periodic orbit. One stretch is 
necessarily close to the periodic orbit, but when the other 
also becomes short and correlated with the central peri- 
odic orbit, it too must follow the periodic orbit closely. 
The orbit 7 then only encounters the periodic orbit once, 
follows it for some number of traversals and then exits 
the system. With no way to swap traversals between the 
different encounters, the partner 7' is then identical to 7 
and included in the diagonal approximation. 

Orbital configurations with periodic orbit encounters 
as described above also occurred in the calculation of 
the covariancc of the reflection coefficients il9<3 yielding a 
term proportional to (l — e~ 2TE / TD ). In fact, by cutting 
the top loop in Figure 0] and moving the ends to the cor- 
rect places we can see we recreate Figure Reversing 
this cutting though, to return to the transmission and 
reflection, we create the second periodic orbit which is 
the top loop in Fig. @] and travels through f 3 , the en- 
counters and £4 in Fig. [3] We will see that this changes 
the orbital configurations compared to the case of the 
variance, changing also the resulting contribution. For 
the covariance of the reflection coefficients it turned out 
to be essential [l9| to consider additionally to the en- 
counter stretches, where both orbits in Fig. [5] are corre- 
lated with the central dashed-dotted periodic orbit, also 
the encounter fringes, where the orbits are correlated 
with themselves (or each other) but where they are no 
longer correlated to the periodic orbit. We marked the 
places where correlations between fringes occurs in Fig. [2] 
by black vertical lines. The duration of the fringes be- 
fore the orbits get correlated to the central periodic orbit 
is denoted by t s and after the orbits leave the central 




FIG. 4: Diagram studied in the calculation of the first quan- 
tum correction to the transmission and reflection in the ab- 
sence of time reversal symmetry. In this example, the or- 
bit traverses the central periodic orbit (dashed-dotted line) 
once during its first encounter and twice during its second 
encounter. We draw the parts of the orbit during the sec- 
ond encounter dashed to distinguish them from the first. The 
partner orbit (not shown) has one traversal of the central 
periodic orbit exchanged between its first and the second en- 
counter with the periodic orbit (i.e. it goes around twice then 
once). The fringes are marked by black vertical lines perpen- 
dicular to the trajectories. 



periodic orbit by t u . 

These fringes are the key to the difference between the 
possible orbital configurations for the covariance of the 
reflection coefficients on the one hand and the transmis- 
sion and reflection on the other hand: In the case of 
the covariance of the reflection coefficients these fringes 
need to have a nonvanishing length, because the two or- 
bits (see Fig. [2J which are correlated during the fringes 
have to end at two different leads where they have to 
be uncorrelated. The orbits away from the central peri- 
odic orbit must be long enough for the chaotic dynamics 
to allow this to happen. When we join one end of the 
dashed and one of the solid orbit in Fig. [5] say to return 
to the transmission (or reflection) as in Fig. [4] then it 
is no longer necessary that the upper periodic orbit (in 
Fig. m thereby created has to be longer than the fringe 
times. These fringes can now start to overlap as depicted 
in Fig. [5j compared to Fig. [3] we let the fringes grow 
till they overlap in the link £3 which itself is of positive 
duration. The stretches of the orbit that connect to the 
leads must though still be longer than the duration of the 
fringes as they must decorrelate from the upper periodic 
orbit to exit the system. 

Now we can explain the effect of the possible orbital 
configurations on the resulting contributions. For this 
we first review some details of calculations for obtaining 
contributions from orbits differing in encounters from 6, 
[T3 |. In order to count the number of orbits we use a 
sum rule based on the classical ergodicity which takes 
the form [l2| 

poo 

\B 7 f = / p(t) (8) 
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FIG. 5: Example of an orbit, that has to be considered in the 
calculation of the transmission, however not in the case of the 
calculation of the conductance fluctuations. The parts of the 
orbit that were changed in comparison to Fig. 0] are shown 
dotted (red). 



with N\ and N2 the number of open transverse chan- 
nels in the left and right lead respectively and p(t) the 
survival probability of an orbit of duration time t. This 
probability decays exponentially for chaotic systems as 
p(t) ~ e~ t/TD for t ->■ oo with t d = T H /N and N = 
N% + N2- The encounters are characterized by the sta- 
ble and unstable coordinate differences s, u in a Poincare 
surface of section inside each encounter between the cen- 
tral periodic orbit and the surrounding orbit. In terms of 
these coordinates the action difference between the two 
trajectories is given by 0, HU 



(9) 



where the sum runs over the different encounters of the 
considered orbit with the central periodic orbit. The 
length of each encounter is obtained by choosing that 
the considered orbit and the central periodic one close 
enough to each other that they can be linearized around 
each other Q 



i(Sj, : Hi) 



hi 



(two connecting the opening to the central periodic orbit 
and one the periodic orbit to itself), r p is the duration of 
the central periodic orbit and t' the time difference be- 
tween the two points (in the different encounters) where 
each of the two encounter stretches reaches a phase space 
difference c with respect to the periodic orbit. The lim- 
its of the time integrals in (fTT]) are determined by the 
fact that the duration of the links, the periodic orbit and 
the encounters have to be positive. This differs from the 
treatment of Fig. [3] in [l3| as they assumed that all five 
links have to have positive duration, but we allow some 
of them to overlap. This automatically includes the other 
cases described at the start of this section as part of a 
continuous deformation of Fig. [3] or Fig. |4] In particular 
we allow ti and £5 to shrink and instead just assume that 
t p is positive, so we therefore use this variable in 

Using these quantities, we obtain from the definition of 
the transmission T (J7]) the following contribution result- 
ing from th e dia grams shown in Fig. 2] and Fig. [SJ which 
we denote T^El 13] 



N1N2 



dt 



d 2 sd 2 uw t {s,u)e^^ As p'(t) 



(10) 



We defined here the modified survival probability in the 
presence of encounters p'(t) taking into account the mod- 
ification mentioned after Eq. (|10j) for encounters within 
the fringes and with periodic orbits: when the encoun- 
ters surround the periodic orbit the parts of the encounter 
stretches traversing a certain point of the periodic orbit 
are so close to each other that they either leave the cavity 
during the first traversal or do not leave at all [lH [l9| 
leading to 



p'(t) — p(i) e ( t ™ c ' 1+t <=° c . 2+ta+ *")/ T£5 = e -(*l+*2+*3+n>)/TD_ 

This expression can be transformed, using and con- 
verting the integral over the full duration of the orbit t 
into one over the link £3, into 



with c of order one. During t enci the survival probability 
is enhanced: either the orbit leaves the system during 
the first stretch or does not leave at all. A density of en- 
counters Wt{s, u) with respect to s, u which characterizes 
the expected number of encounters the orbit has with the 
central periodic orbit can then be obtained, in the case 
here of two encounters with the periodic orbit, as 13] 



w t (s,u) 
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enc,l^enc,2 




n/* 



where the phase space volume of the system under con- 
sideration is denoted by Q and ti denotes the duration 
of two of the three links away from the periodic orbit 
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(14) 



i=l 



where we also used the explicit form of the survival prob- 
ability and the action difference AS*. To understand that 
the expression in (|14[) yields zero, we perform the inte- 



grals with respect to Sj,itj, like in [li 



d Sl du t e (t/h)s * u > — 



1/Xi 



2 

C Xj 



l 



= 4c 2 / dxi I doi cos i,i 

1 J 1 V " / 0°i*enc,i 

4c 2 A / cix,- cos 



(15) 



with the substitution it; = cjoi and = cXiUi. The 
integral in the last line in (|T5|) rapidly oscillates as a 
function of energy in the limit H — > and thus yields no 
contribution due to the energy average in (fl"4")l . 

We thus obtain that there are no quantum corrections 
(at least to this order) to the transmission when time 
reversal symmetry is absent 



0. 



(16) 



and a similar calculation shows that this also holds for 
the reflection R. Coherent backscattering, i.e. having 
encounters at the opening that additionally have to be 
taken into account for reflection, is also not possible. 
First this requires the encounter to be traversed in op- 
posite direction on both traversals, which can only occur 
with time reversal symmetry. Second, even with time 
reversal symmetry, when the trajectory returns to the 
encounter the second time it would necessarily escape 
the systems, and not be able to complete the rest of the 
semiclassical diagram. 

To summarize, we saw in this section how, despite their 
close similarities, the two different orbital configurations 
appearing in the case of the covariance of the reflection on 
one hand and the transmission and reflection coefficients 
on the other lead to two different results: in the case of 
the covariance of the reflection to a term proportional 
to (l — e~ 2rE / TD ) , in the case of the transmission and 
reflection coefficients to zero contribution. 



B. With time reversal symmetry 

We now turn to the calculations in the case with time 
reversal symmetry. In this case we also have to consider 
diagrams where the encounters are traversed in different 
directions by the orbit. As their contributions are quite 
different we will study them individually. We start with 
two independent encounters with no central periodic or- 
bit involved, referred to as two 2-encounters, shown in 
Fig. [5] We first cut the periodic orbit during one of the 
middle links and refer to the corresponding contribution 
as T^t. In this case the contribution of the two s,u- 
integrals for the two different encounters factorizes and 
can be evaluated for each encounter separately, as was 
done in [l9| for the case of the reflection covariance (ob- 
tained by cutting both the leftmost and rightmost links 
in the periodic orbit in Fig. [5]). Each encounter provides 




FIG. 6: Periodic orbit with two independent 2-encounters. 
The different positions, where it can be cut to obtain an 
open orbit contributing to the transmission are indicated by 
(red) perpendicular lines, the position of the (blue) encounter 
stretches are indicated by a box. 



a factor ~Ne t e/ t d ^ the five links factors N 1 and the 
leads the factor NiN 2 so that we obtain for the contri- 
bution I®\ 



Tit = ^1^2 c -2^/r n 

(iVi+7V 2 ) 3 



(17) 



The corresponding contribution to the reflection iiP 1 is 
obtained by multiplying I® 1 by N 1 /]\J 2 to take into ac- 
count that the orbit leaves through the lead 1 instead of 
lead 2. When cutting the left link of the periodic orbit 
in Fig. |6l whose contribution we denote T^, we obtain 
for the transmission the same result as I® 1 . However 
for the reflection in this case it is also possible to obtain 
a coherent backscattering contribution by shrinking the 
length of both links on the left in Fig. [5] to zero (or we 
cut the diagram in Fig. [6] at the leftmost encounter and 
move this to the lead). Also in this case the encounter 
integrals for the two encounters factorize, yielding 



TV 2 



(iVa + N 2 



-2te/td _ 



(JVi 



7C 



-2te/td 



N 2 ) 

(18) 

where the first term is the same as fP* 1 and the second 
comes from the coherent backscattering. 

Next we consider, as for the case of no time rever- 
sal symmetry, the situation of two 2-encounters near a 
periodic orbit. The configuration where the encounter 
stretches are parallel (in the same direction) was treated 
in the last subsection so, as we have now the freedom 
to traverse the two encounter stretches in opposite direc- 
tions, we now turn to configurations where some of the 
stretches are antiparallel to each other. Starting with the 
periodic orbit configuration in Fig. [7] (a 3-encounter in 
[g]) there are three possible places to cut this orbit open 
as shown by the red lines perpendicular to the orbit. By 
opening the parts not enclosing the central periodic orbit, 
we obtain a configuration shown in Fig. [SJ 

Unlike the case without time reversal symmetry, we 
can see that some of the fringes must have non- vanishing 
length like considered in [19J: the two fringes marked by 
dotted boxes in Fig. [5] on the right hand side of the en- 
counter cannot have vanishing length, because as long as 
the two parts are correlated, the corresponding loop they 
form cannot close. The two fringes in Fig. [5] on the left 
hand side of the encounter can only vanish in the case of 
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FIG. 7: A periodic orbit encounter only existing in the case of 
time reversal symmetry. The central periodic orbit is drawn 
with a dashed-dotted line, the position of the (blue) encounter 
stretches are marked by a box. The (red) lines perpendicular 
to the orbit mark the places where it can be cut open. 




FIG. 8: A periodic orbit encounter only existing in the case of 
time reversal symmetry. The position of the (blue) encounter 
stretches are marked by a box, the position of the fringes by 
small dotted boxes. 



coherent backscattering, i.e. if the orbit starts and ends in 
a correlated manner in the same lead. Note that in the 
left fringe (defined where stretches are correlated with 
each other away from the central periodic orbit) we only 
have the two stretches which connect to the leads and 
that the remaining encounter stretch in Fig. [5] which fol- 
lows the central periodic orbit has already decorrelated 
from the others so it does not also need to escape. To 
evaluate the contribution we first need to determine the 
values of the prefactors a, b, d in the exponential in J in 
the Appendix in front of i e nc,i +ienc,2, t a +t u and t p , re- 
spectively. As the survival probability along the periodic 
orbit depends only on r p and not on i nc,i> tenc,2) we ob- 
tain a = and d — 1/t£>. During the fringes we have two 
correlated stretches with the survival probability deter- 
mined by one of them, thus yielding b = — 1/td. When 
multiplying the resulting contribution for J by the fac- 
tors resulting from the links and the channel factors due 
to the leads we obtain the contribution 7^ originating 




FIG. 9: A 3-encounter involving no periodic orbits. The 
(blue) encounter stretches are again marked by boxes, the 
places where fringe correlations can occur (are marked red) 
indicated by smaller dotted boxes. 



from Fig. [5] to the transmission 



NlNo 



^—2t e /t d 



2(N 1 + N 2 y 
and to the reflection 



(19) 



2(N 1 +N 2 y 
2 (JVi + N 2 f 



1 _ g-2r E /T D 



I _ c -2r E /r D 



(20) 



The last case, depicted in Fig. [9l is obtained by open- 
ing along the central periodic orbit in Fig. [7] In general 
any two of the three stretches on either side of the en- 
counter could remain correlated in the fringes away from 
the main encounter where all three stretches are close 
and correlated. The duration of the fringes, i.e. here in 
general the orbital parts where only two of the three en- 
counter stretches are correlated, is denoted before and 
after where all three orbits are correlated by t s and t u , 
respectively as in [21] and in the Appendix. On each 
side, fringe correlations become important if the two or- 
bital parts containing the fringes are connected to each 
other, referred to as case A, but not if one orbital part of 
them is connected to the opening, referred to as case B. 
The reason why we have to take into account fringe cor- 
relations in case A is that namely the loop cannot close 
as long as the two parts of the orbit are still correlated. 
In case B the part of the orbit connected to the open- 
ing still has to be longer than the fringes so that when 
it escapes it does not force the rest of the orbit to also 
escape. However the other fringe which lies on the orbit 
that is not connected to the opening in Fig. |H] has no 
length restriction: if the length of that fringe tends to 
zero, the orbital part connected to the opening will just 
follow the first one for the time t s or t u . The latter part 
then also contains the survival probability contribution 
due to the fringes. The fact, that only the part of the 
orbit connected to the opening has a length restriction 
due to the fringes together with the enhancement of the 
survival probability during the fringe parts lets, as al- 
ready in (|14[) . the t s and the t u drop from the resulting 
expressions for the contribution 2^ from these diagrams 
in case B. 
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With these remarks in mind we evaluate the contribu- 
tion in Fig. IHlby making use of the results obtained in 
pll ] , that we review in the Appendix as contributions K i 
and K 2 . The overall contribution K is split into two parts 
K\ and K 2 : K\ contains the contribution resulting from 
the 3-encounter without fringes and K 2 the contribution 
resulting from the difference between the 3-encounter 
with fringes and a 3-encounter without fringes. In all the 
cases considered here we include the first part K\ where 
as for the survival probability during encounters we only 
need to count one encounter stretch we get / = — 1/td 
in K\ in the Appendix. To obtain the contribution K 2 in 
this case we first note that it was shown in [2l| to be suf- 
ficient to only consider certain encounter diagrams: only 
one stretch contains two fringes, the other two fringes lie 
on the other stretches in a certain way. Furthermore by 
setting gx = or g 2 = 0, implying that only t s or t u is 
nonzero we obtain zero contribution, see Eq. (|55|) . Thus 
three different possibilities remain: one from each of the 
three stretches containing two fringes, one belonging to 
case A and two to case B. As already explained there 
is no t a , ^-dependence in case B and thus in this case 
no contribution to K 2 . In case A we obtain 1/3 of the 
contribution K 2 in Eq. (|54|) . we set / = g = — 1/td to 
take into account that only one stretch of the encounter 
is taken into account in the survival probability. For the 
overall contribution we therefore have 



mN 2 



(jvi + N 2 y 



-2te/td 



(21) 



The same contribution times the factor Ni/N 2 is ob- 
tained for the reflection. 

To calculate the overall quantum correction to the 
transmission at the considered order T 2n a in the case of 
time reversal symmetry we sum twice (to account for dia- 
grams related by symmetry) the contributions from (|17[) 
and (fT9l) and the contribution from ([21]) yielding 



T 2 nd = 



NN 2 



(Ni + N 2 



(22) 



Note that this quantum correction is independent of the 
Ehrenfest time. This also holds for the corresponding 
contribution to the reflection i?2nd, which we obtain here 
by adding the contribution from (|18|) to twice the contri- 
bution from ([20} and to the related contributions from 
and fl2U) multiplied by N x /N 2 



-ffend = 



(m + N 2 y (Ni + N 2 y 



Current conservation 



(23) 



Having calculated all contributions to the transmission 
and reflection we now want to check if current conserva- 
tion is fulfilled, i.e. if the transmission and the reflection 
calculated for one lead add up to the number of open 



channels in that lead. As without time reversal symme- 
try there are no contributions at the order 1/N consid- 
ered here, current conservation, already fulfilled at the 
diagonal level, is thus not violated. In the case of time 
reversal symmetry the contributions to T and R at the 
considered order are given in (|2"2"j) and (|2"3"]l and sum to 
zero. Current conservation is thus again fulfilled. We 
want to emphasize here that correlations between en- 
counter fringes, first treated in 19], were important to 
obtain this result: forgetting for a moment the effect of 
fringe correlations, the contribution (f2~Tj) would possess 
the Ehrenfest-time dependence cT TE l TD and the contri- 
butions (119I20P would be zero leading to a non current 
conserving result for T 2n a and i?2nd- 



III. FREQUENCY DEPENDENCE 

In this section we want to generalize the results ob- 
tained for dc-transport to the ac-case (26|, i.e. we want 
to consider 

T(lu) = (Tr(t(E + tvjj/2)t\E - hu/2))) (24) 

and a correspondingly defined R(uj). As the calculation 
leading to this generalization is straightforward we only 
briefly explain the difference to the calculation before and 
then show the results. In general adding a frequency de- 
pendence means including into the formulas in section 
II a factor e lut with the overall duration t of the or- 
bit. In the case of no time reversal symmetry we get 
in terms of the notation of eq. (TT4"} an additional fac- 
tor e^h+'^^+^-^ntie""' 1 . To include this factor 
when performing the s, u- integrals we take a, b, d for 
lj = from the last section and include the w-dependent 
exponential factor given in the last sentence to obtain 
a = iuj, b = and d = 1/t£> — ilj. Inserting this in J in 
the Appendix and taking into account the factors from 
the links and the leads we obtain 



-NxN 2 



(Nl + N 2 f (1 - lLOT D f 



In the orthogonal case, including a frequency de- 
pendence into eq. (|17|) adds an additional factor 

p icj(2t e „ c , 1 +2t e „ c , 2 ) 



it; 



yielding finally 



= 2»( W ) = ^|iP» 

7Vi7V 2 (1 _ 2iLOT D f 



(A!+7V 2 ) 3 (l-iU!T D f 



-2te/td+4iutb 



(26) 

The first term of eq. (fT5]) is modified in the same way 
as the expression in eq. (|26l) while for the second we 
have three links instead of five, reducing the power of 
(1 — iojTjj) in the denominator by two, and an additional 
integral over the duration of the encounter reducing the 
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power of (1 — 2iuTu) by one compared to the first term. 
We thus obtain 



(1 - 2iujT D ) e _ 2TE/TD+4iuTE 

[Nl+N 2 f {l~iiJT D f 



Ni 



(1 — 2iujtd) 



{Ni + N 2 f (1 - ilutd) 3 



-2T E /To+ii^T E 



(27) 



In equation p9j) the additional factor 
e iu(T p +t cnc , 1+ t cnc , 2 +2t 3 +2t u ) lf i=1 e^ occurs, the equa- 
tion is thus replaced by 



7%) = 



f g2iii)TB _ g— 2t e /t d +Hu)t e 



2 (Ni + NaT 
(1 — 2iujTE>y 
(1 - iu)T D f (1 - ilotd) 



2lu 2 t d 



(28) 



The latter equation can be obtained from J in the Ap- 
pendix by setting a = iu, b = —X/td + 2ico and d = 
1/td — iw, again considering the additional terms from 
the w-dependent exponentials. The additional frequency 
in the first term in eq. (|2"0")) has the same effect as in eq. 
(fT9"| . in the second term we again have one instead of 
three link times i, and an additional integral over t s or 



N 2 



2{N 1 +N 2 ) 
(1 — 2iu)TnY 
(1 - iuitd) 

2(N 1 +N 2 f 

(1 — 2iojtd) 
(1 - iu)T D ) 



(1 - iuitd)" 



' g2itJTE _ 2t e /t d +Hut e 



(29) 



In case of eq. (f2T|) we get by taking the corresponding 
contribution of the encounter again from the Appendix 
with / = — 1/td + 3ia> and g = — 1/td + 2iw in i^i and 
i^2 since £ e nc is traversed three times and the fringes two 
times 



oj = -- 



A^iA 2 



[Ni + N 2 y 



(i - iu>T D y 



(1 — 2iujT£,) 2 

(1 - iU!T D ) 5 
-T E /r D +3iuiTE 



3 — 2te/td +Huite 



(30) 



and a corresponding contribution for the reflection. 

After obtaining these results it is now possible to check 
if they fulfill the relation 



-4- \ S ( E + S f [E - hLo)} 



= 



(31) 



with the scattering matrix at the energy E, S(E) con- 
taining the reflection and transmission subblocks r, t for 
the incoming wave in the lead 1 and r' , t' for the incoming 
wave in the lead 2, respectively 



S[E) = 



r[E) t'[E) 
t[E) r'[E) 



(32) 



Before we only considered the reflection and transmission 
for an incoming wave in the lead 1, i.e. only the corre- 
lators of elements of r[E) and t[E). The corresponding 
results for the correlators of r'(E) and t'[E) are obtained 
by swapping Ni and N 2 . 

In order to see why relation (f5T|) is fulfilled we 
rewrite it in terms of the Wigner time delay [27| . 
measuring the additional time spend in the scatter- 
ing process compared to the free motion, = 



±Tr [S [E + hcu) S{E- Hu))] 



-T W = 0. 



Equation (pTTj) is then 



(33) 



That this relation has to hold can be obtained by compar- 
ing the two equivalent representations of the Wigner time 
delay discussed in [27| : their semiclassical equivalence is 
discussed in 22j. The first representation in terms of 



the density of states involves a single sum over trapped 
periodic orbits, the second representation in terms of 
transmission coefficients involves a double sum over lead- 
connecting paths. As the first representation yields, after 
taking an energy average, an Ehrenfest time independent 
result - we cannot identify any Ehrenfest time dependent 
contributions in a single sum over periodic orbits - tw 
has to be Ehrenfest time independent. 

In terms of the subblocks of S(E) introduced in eq. 
(1521 . Tw can be expressed as 



duj 



\[T[u) + R(cj) + I» + #H]Uo (34) 



with the primes again denoting that the incoming wave is 
in lead 2 instead of lead 1. We start our further analysis 
of the first two terms: In order to check if our results for 
T(lj) and R(uj) given above eq. (13T1) fulfill relation (|33l) . 
we first consider the sum of the contributions to tw which 
decrease with increasing Ehrenfest time. The contribu- 
tion proportional to e~ TE / TD is obtained by considering 
the corresponding term in (|30[) yielding 



d_ 



Ni 



[Ni+N 2 f {1-ilutd) 5 



=0 

(35) 

For calculating the contribution proportional to e 2te I Td 
we sum the corresponding terms from (|2f)|) . (|27l) and ([50)1 



d 

did 



N x uj 2 T 2 D (l-2iwT D ) 
(Ni+N 2 ) 2 [l-iujT D f C 



2te I 'td+Huite 



oj=0 



(36) 
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In the case of the contributions increasing or oscillating 
with increasing Ehrenfest time we obtain from ([25)1 and 
(EH) 



2^2 



Ul T 



(JVi + N 2 f (1 - ilotd) 5 



[2e : 



— (1 — 2iLJT]j) ^( 

and from (1251) 



2i<jJTE 



-2te /td+Muite 



)]} (37) 



N 2 f (1 - iLOT D f 



(38) 



U1 = 



which is the only contribution also existing in the absence 
of time reversal symmetry. 

The results in (|35H38p fulfill eq. (13"TT) . because all are 
proportional to oj 2 and thus are equal to zero after dif- 
ferentiating them with respect to u> and setting uj = 0. 

The results obtained from the two second terms in eq. 
(fM)) differ from the first ones by a factor N 2 /N\ and thus 
also yield zero contribution to tw- 



IV. SPECTRAL FORM FACTOR 

In this section we want to apply our knowledge about 
the orbital configurations which contribute to the con- 
ductance to calculate the first off-diagonal quantum cor- 
rection to the spectral form factor. We first want to 
briefly review the contributions calculated in 23] so we 
will use almost the same notation as there and for further 
details we refer the reader to that paper. For systems 
with time reversal symmetry the first correction derives 
from the orbit pair depicted in Fig. Q] whose Ehrenfest 
time dependence is simply e~ TB ^ TD . For systems with- 
out time reversal symmetry however we have diagrams 
starting like in Fig. EH but with the orbits in the leads 
connected together so that t\ and t 2 join to a single link. 
Note that we can then identify four periodic orbits in 
the picture, one central orbit, one through t% as before 
and two through the newly joined links and t% and £5, 
respectively. From there we can allow the encounters 
to overlap to create a 3-encounter and then finally to 
wind around the central periodic orbit, as described at 
the start of section III Al As we saw for the transmission 
there are further possibilities compared to the covariance 
of the reflection (or we can relax more restrictions) and 
likewise here there are additional contributions. We will 
see how they lead semiclassically to the field-theoretical 
result for the first off-diagonal quantum correction to the 
spectral form factor, but first we recall the results of the 
diagrams covered in [23[. 

In [23j ] the contribution of two independent 2- 
encounters (c.f. Fig. [3]), denoted by SK 2 i,(uj), is given by 



SK 2b (uj) 



1 



Q2 e 4iwT E 



2-kKT^ Ouj 2 oj 2 



(39) 




1Z 



V 



FIG. 10: A diagram accounted for in the contribution 
SK2d(u) to the spectral form factor. A central dashed-dotted 
periodic orbit is encountered two times. Fringe correlations 
are marked by black vertical lines. For the partner (not 
shown) one traversal of the central periodic orbit is exchanged 
between the first and the second encounter. 



the contribution of one 3-encounter, 5K 2c (uj), obtained 
by allowing the encounter stretches to overlap along one 
enclosed periodic orbit in Fig. [3] at one end, is obtained 
using [21| to be 



SK 2c (cj) = 



1 



d 2 i 



2-kHT% 8uj 2 oj 2 



3e 



3iCJTe 



4e 



4i0JTE 



) ■ (40) 



A further diagram results from encounter overlap along 
one enclosed periodic orbit at both ends, see Fig.HUl The 
overall contribution / containing the contribution from 
the latter diagram and the contributions SK 2 b(uj) and 
8K 2c {uj) is obtained by considering J in the appendix 
with a — iuj, b — 2iui and d = —iuj, because the orbit 
is assumed to be longer than t e nc,i + icnc,2 + 2t s + 2t u , 
and multiplying it by factors resulting from the links not 
surrounding the central periodic orbit. A technical com- 
plication is that this diagram contains three copies of the 
contribution with a 3-encounter and 4 copies of the con- 
tribution with two 2-encounters. Naturally we only want 
to include one copy later so we subtract them all here. 
All told, the contribution resulting from theperiodic or- 
bit encounters, 8K 2 d (w) was calculated in [23j to be 



SK 2d (uj) = I-4SK 2b {uj)-35K 2c (uj) 
1 d 2 1 



2nhT 2 duS 1 uj 2 



3e 



2icJTE 



9e 



+6c 



(41) 



However, as already explained above, this is only the 
complete set of contributions in the case when the two or- 
bits approaching and leaving the periodic orbit are open 
like in the case of the reflection covariance as in Fig. [5J 
otherwise the corresponding orbital parts can get shorter 
than the duration of the fringes, like moving from Fig. 0] 
to Fig. [S] To take into account this additional configura- 
tion we replace SK 2c i(uj) by another contribution denoted 
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by SK2e (w)- Up to now we only allowed the orbital parts 
decorrelated from the central periodic orbit, to be longer 
than 2t s + 2t u . However the contribution considered now 
results from an orbital configuration where the two other 
links, those decorrelated from the central periodic orbit 
in Fig. [lOj get shorter than t s + t u each. The bottom 
and the top loop in Fig. [TU] outside of the encounter with 
the central periodic orbit must again have positive length 
but do not necessarily need to be longer than the fringes. 
In order to calculate this contribution we first consider 
J in the Appendix with a = ilo, b = and d = —iuj, 
because the orbit has a minimal length of i e nc,i + t cnCi 2 
as in Eq. (|25|) , together with the factors resulting from 
the links not surrounding the central periodic orbit and 
subtract from this contribution like in Eq. (|41|) 5K2b(u) 
and SK^ciuS) with the right multiplicity factors. The cor- 
responding contribution denoted by 5K2 e ,i{oj) is 



8K 2e ,\{w) = 



1 



d 2 1 



2irM% du 2 lo 2 



4<5K 2b H-3<5K 2 c(w) 



d 2 1 



2ttHT 2 8uj 2 lo 2 



9e 



8e 



(42) 



This procedure however counts some configurations con- 
taining a surrounded periodic orbit twice: shrinking 
in Fig. [TUJ the length of the upper periodic orbit to 
zero we again obtain a contribution containing one sur- 
rounded periodic orbit. A configuration containing one 
surrounded periodic orbit was however already taken into 
account in 8K2 e ,i{w) when shrinking the length of the 
central periodic orbit to zero. We thus subtract the lat- 
ter contribution. This contribution is calculated by again 
making use of K\ and K2 in the Appendix: we therefore 
consider a 3-encounter, i.e. / = 3iw, with fringes with 
duration between t s + t u and 2t s + 2t u . We thus con- 
sider once the prefactor g — 2iuj and once g — iuj in front 
of t s + t u and take the difference of the two results ob- 
tained for K in the Appendix yielding the contribution 
SK 2 e,2( UJ ) given by 



6K 2e ,2( UJ ) 



1 



d 2 1 



2irhT 2 duj 2 lo 2 



4 



)] 



(43) 



Adding the two contributions to <5ir 2e (cj) we obtain 
1 d 2 4 



2nhT 2 duj 2 lo 2 



( e 4i ^ - e 3iwrE ). (44) 



Summing now the quantum corrections in the absence 
of time reversal symmetry given in (|39|) . (j40|) and (j44|) 
we obtain for the overall quantum correction at the con- 
sidered order 5K(lu) 



6K(w) 



1 d 2 (e 4i " TE — e 3iwrE ) 
2irhT 2 9^2 -ji 



(45) 



This yields then after the Fourier transform for the cor- 
responding correction to the spectral form factor SK(t) 

T 2 

5K{t) = -— [9 {tT h - 3t e ) - 6 (tT h - 4t e )} (46) 

with <d(x) = L dx'9(x') = xO{x) with the Heaviside 
theta function 9(x). Expression (l46l) was also obtained 
in [23] by field-theoretical methods. 

As already noted for the conductance, we also want to 
emphasize here that these results for the spectral auto- 
correlation function could not be obtained without con- 
sidering fringes: not doing so we would only get the con- 
tribution (|39[) with a multiplicity factor four along with 
the contribution from a 3-encounter with 3 equally long 
encounter stretches, given by 



q2 e 3iuT E 



(47) 



with a multiplicity factor three and the overall contribu- 
tion resulting from all possible encounter configurations 
given by 



1 



2nhT 2 du; 2 uj 2 



(48) 



with a multiplicity factor one. As one can easily see it 
is not possible to obtain the field-theoretical result from 
just these semiclassical contributions. 



V. CONCLUSIONS 

In this paper we have shown how to calculate the (1/iV) 
quantum correction to the transmission and reflection for 
systems both with and without time reversal symmetry. 
Starting with dc-transport, we obtained at the consid- 
ered order that the transmission as well as the reflection 
are zero in the case of no time reversal symmetry. In 
the presence of time reversal symmetry the overall con- 
tributions to the transmission (j2"2"j) and the reflection (f2"3")l 
are independent of the Ehrenfest time and fulfill current 
conservation. This simply means that the quantum cor- 
rections to the transmission and the reflection add up to 
zero. We extended this analysis then to the ac-transport 
by including a finite energy difference frw between the two 
scattering matrix elements. For the Wigner time delay 
we saw that the results led to zero extra contribution and 
importantly that there is no Ehrenfest time dependence 
consistent with the two complementary semiclassical rep- 
resentations of the time delay. 

For the transmission and reflection the key step is that 
we can relax one of the restrictions compared to the cal- 
culation of the reflection covariance in [19j . Namely, be- 
cause of the slightly different topology formed by rejoin- 
ing some of the links (previously cut to get to the reflec- 
tion covariance) the fringes are allowed to overlap and 
surround the second periodic orbit formed. For closed 
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systems without time reversal symmetry we then rejoin 
more links and create a third (and fourth) periodic orbit. 
Also this relaxes a restriction and leads to a modifica- 
tion of the results obtained for the spectral form factor 
in 23] . By including all possibilities we showed that our 
semiclassical result agrees with the field-theoretical pre- 
diction from [24j . and hence provides some justification 
for the phenomenological treatment of effective RMT. 

These results are a first step towards the semiclassical 
calculation of transport and spectral properties including 
the Ehrenfest time dependence to arbitrary high orders. 
In the case of transport we found that the second order 
quantum correction is independent of the Ehrenfest time, 
whereas the first quantum correction turned out to be 
proportional to q~ te / td [l7j : it seems thus not yet pos- 
sible to identify a general pattern behind the Ehrenfest 
time dependence of the total contribution to this trans- 
port quantity at different orders. For the related case of 
the spectral form factor the end result seems to possess a 
simple structure: the contribution from each discrete di- 
agram (i.e. those considered without the Ehrenfest time 
dependence in Q) to the spectral form factor K{t) is 
given by a function depending on t — tite where n is 
the number of encounter stretches of the underlying dia- 
gram. This is the structure predicted by effective RMT 
and the next step here would be to see if this indeed 
holds semiclassically to arbitrarily high order. In princi- 
ple this would require a general treatment of arbitrarily 
sized periodic orbit encounters and their fringes and we 
have the main ingredients [19l |21| but lack the diagram- 
matic rules, that are available when the Ehrenfest time 
vanishes [lj, [l5| , for treating all possible correlations oc- 
curring mainly during fringes. We wonder if it would be 
possible to partition the semiclassical diagrams in an ef- 
ficient way so that any such structure becomes clear, as 
was the case for the related problem of the correlator of 
2n scattering matrices (28j . 
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VII. APPENDIX 

In this appendix we want to calculate the integrals over 
the stable and unstable coordinates s and u occurring 
for encounters with periodic orbits for general prefactors 



a, b, d in front of i e nc,i +Wc,2, t s +t u and Tp, respectively. 
We therefore consider the expression J defined as 



dr p I d sd u I dt' 



1 



I ! J —c JO ^ 2 ^cnc,l^cnc,2 

x exp (a (i cnc ,i + t enc<2 ) +b(t s + t u ) - dr p ) 
X exp I - ^ s t Ui \ (49) 



with the duration of the fringes before and after the con- 
sidered orbit approaches the central periodic orbit, t s and 
t u , respectively. The s,u- and the i'-integrals in J are 
evaluated in [19| yielding 



JO J-H 



oo L2 -drp 

d TP (e 2aTE - e 2bTE ) 



Performing the rp-integral finally yields 



(50) 



a b 



e 2ar E _ e 2br E j ^ 



H 



2Tf I d (a - b) 



which is frequently used in the main part of this article. 

Furthermore we want to consider a 3-encounter and 
calculate with arbitrary prefactors / and g in front of 
t onc and t s + t u respectively 

K = (J d 2 sd 2 u-^ — exp (ft cnc + g x t s + g 2 t u ) 



xcx p( IE' 



(52) 



The first contribution K\ is obtained by setting g = 0, it 
yields |l| 



= — e /TE 



T 



(53) 



H 



The second contribution Ki given by K — K\ is given for 
.91 = 92 = g by [HI 



^2 = 3- 



'(2g-f)T 2 
For gi ^ <?2 we have instead [H 



2gr E _ „/te^ 



K 2 = 3 



5i 92 



(.91 + .92 - f)T 2 



(V 91 



+ 32)l"E _ p /l"E 



(54) 



(55) 
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